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Abstract 

We study the stabihty of network communication after removal of g = 1 — p links under the 
assumption that communication is effective only if the shortest path between nodes i and j after 
removal is shorter than alij{a > 1) where £ij is the shortest path before removal. For a large class 
of networks, we find a new percolation transition at pc = (kq — 1)(^~")/", where Kq = {k'^)/{k) 
and k is the node degree. Below p^, only a fraction A^"^ of the network nodes can communicate, 
where 5 = a(l — | logp|/log (kq ~ 1)) < li while above Pc, order nodes can communicate within 
the limited path length aiij. Our analytical results are supported by simulations on Erdos-Renyi 
and scale-free network models. We expect our results to influence the design of networks, routing 
algorithms, and immunization strategies, where short paths are most relevant. 
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The study of complex networks has emerged as an important tool to better understand 
many social, technological, and biological real-world systems ranging from communication 
networks like the Internet, to cellular networks In many cases, networks are the medium 
through which information is transported, i.e., in social networks the propagation of epi- 



demies, rumors, etc. and in the Internet the propagation of data packets 

An important question regarding networks is their stability, i.e, under what conditions 
the network breaks down llO|. In communications, a network breakdown means 

information cannot be transmitted to most nodes, and in epidemiology, that an epidemic 
has stopped. 

rn 

The main approach for studying network stability is percolation theory |11]. In perco- 
lation, a fraction q = 1 — p of the network nodes (or links) are removed until a critical 
value pc is reached. For p < Pc the network collapses into small clusters, while for p > pc, a 
spanning cluster of order nodes appears @, ol, H, Q, Q. However, even though in the 
original network the nodes are connected through short paths, near pc the paths become 
very long. For instance, in the original Erdos-Renyi network the typical distance between 
nodes is of order logA^ [l^ compared to order A^^^^ near the percolation threshold jl^ . 
These long distances may have a significant influence on network function. For example, in 
communication, long paths are usually inefficient, and in epidemics, disease spreading often 
decays in time due to mutations or natural immunization, so for long paths the epidemic 
may die out before the network collapses. In these cases the interesting question is some- 
times, not when does the network break down, but when the network connectivity becomes 
inefficient. 

To answer this question, we propose a new percolation model which we call limited path 
percolation (LPP). In this model, after removing a fraction g = 1 — p of the network nodes, 
any two of these nodes, say i and j are considered connected only if the shortest path 
between them is shorter than aiij (a > 1), where iij is the shortest path before removal. We 
then ask, given our new limited path constrains, what is the value p at which a spanning 
cluster appears. We find a new phase transition, which depends on a, at pc = Pc{o)i where 
Pc < < 1. For Pc < p < Pc, the LPP spanning cluster is only a zero fraction (fractal) of 
the network, which scales as A^'' (5 < 1). For p > pc the LPP spanning cluster is of order A^. 

For simplicity, we start our analysis with Erdos-Renyi (ER) networks and then argue 
that the theory is also valid in general for random networks. We begin with random removal 
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but extend our considerations to targeted removal on highly connected nodes, and find that 
similar phenomena appears. We support our theory with simulations. 

Erdos-Renyi networks are random networks consisting of N nodes connected 

with probability (p and disconnected with probability 1 — 0. The degree distribution ^{k) 
is Poisson with the form $(A;) = {k)'^e'^^'^^ /k\^ where fc, the degree, is the number of links 
attached to a node, and {k) = Yl^=i k^{k) is the average degree of the network. The typical 
distance between nodes is \ogN/\og{k). 

Next we evaluate Sa, the size of the spanning cluster under LPP. After the removal of 
fraction q of the links, the spanning cluster can be considered tree-like since, up to order N, 
loops are negligible |8|. Thus, Sa can be approximated by 

Sa ~ c(p)[p(fc)]"J^ = c{p)N\ 6 = all- < 1 (Erdos-Renyi) (1) 

V log \^/ / I — I 

where p{k) is the average degree after removal, c{p) = Cop{k)/{p{k) — 1) 15|, and 
alogN/ log (k) is the new tree depth imposed by the limited path length restriction. The 
exponent 6 = S{a,p, (k)) is an increasing function of a, i.e., for larger values of a longer 
paths are valid and therefore more nodes are included in the spanning cluster, leading to a 
higher value of 6. The exponent 6 is bounded below by zero and above by 1, since N is the 
maximum number of nodes available. Setting 6 = 1 and solving for p in Eq. ([1]) we obtain 
the transition threshold 

Pcia) = {k)^ (Erdos-Renyi). (2) 

Figure [1] presents the phase diagram for LPP. For Pc ^ P ^ Pc{0') the spanning cluster 
is a fractal of size A^'' and 6 continuously increases with p. For p > Pc{a), a spanning 
cluster of order exists with path lengths i'^j < aiij. Using the function 1 — Pc{a) we are 
able to calculate for a given value of a, the percentage of links that can be removed before 
the network is no longer connected with effective paths, i.e., shorter than aiij. Note that 
for a —>■ oo, when no path length restriction is imposed, we recover the usual percolation 
threshold Pc{a oo) = pc = Jldl- Equations ([I]) and ([2]) are supported by the 



simulations presented in Fig. M^a) 16|, ll8| . For a summary of the various equations in the 
article, see Table [H 

Our results for the different regimes of Sa can be summarized by the scaling relation for 

p> Pc 

Sa-c{p)N'f(^-^^^ (Erdos-Renyi), (3) 



where P^c is the probabihty of an arbitrary node to belong to the usual percolation spanning 



cluster [llj. The function f{x) scales as x when x <^ 1 and approaches a constant as a; ^ 1. 
In Fig. El^a), we present simulation results for several a and p values for ER networks, 
supporting the scaling form of Eq. 

The theory for LPP can be extended to all random networks with typical distance between 
nodes of order logA^ by substituting (k) with the generalized form (k — 1), known as the 
branching factor, defined by k — 1 = — 1 Isl. Replacing {k) with {k — 1) in Eq. ([1]) we 
obtain the general equation for the spanning cluster size 

Sa cip)iK - ly^^^ = N', 5^aM^^^^ (4) 

\og{Ko - 1) 

where Kg — 1 is the branching factor of the original network and k, — 1 the branching factor 
after removal, which depends on p. When a random fraction of the network is removed, 
K — 1 = p{k,o — 1) [s^. Note that for the specific case of ER networks, k, — 1 = p{k) and 
Ko — ^ = (k) , reducing Eq. (j4]) to Eq. ([1]). In the general case of random networks, the LPP 
transition is found by imposing 6 = 1, which yields 

Pc{a) = {ko - '^)^ ■ (5) 

The scaling form for Sa is the same as Eq. with 6 taken from Eq. (jl]). 

Our general theory for LPP can be illustrated on scale-free (SF) networks. Scale-free 
networks have generated much interest due to their relation to many real-world networks. 



such as the 
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nternet, WWW, social networks, cellular networks, and world-airline network 



22|. Scale- free networks are characterized by a power-law degree distribution 
~ fc"'*' {m < k < K), where K = mN^^^'^^^^ [§1. The power-law distribution allows 
a network to have a few nodes with a large number of links ( "hubs" ) which usually play a 
critical role in network function. Calculating k for SF networks one obtains 

For A > 3, Eq. (jl]) is valid and thus LPP is similar to ER networks, except that it depends 
on K — 1 instead of {k). The phase diagram of SF networks is shown in Fig. [U^b). The results 
of the simulations supporting the theoretical value of 5, Eq. (jl]), are shown in Fig. [2]^b), and 
for the scaling form of Sa are presented in Fig. [31(b). 



For 2 < A < 3 the typical network length scales as ^ = 21oglog A^/| log(A — 2)| 23|, 



24l |. For this regime, our scaling approach to calculate Sa is no longer valid since the tree 



approximation breaks down. However, the LPP transition still exists when aiij = •, where 
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i[j is the distance after removal, with typical value i' = 2 log log Pqo^/ I log (A — 2) 
Solving aiij = i'-j for — > oo, we obtain 

« = ^ = ^''f^'f-f - 1 (Scale-free 2 < A < 3). (7) 
iij log log AT ^ ' ^ ^ 

This implies that Pc and thus, for any finite p, Sa is always of order A^. The results of 
the simulations presented in Fig. [2](c) support our prediction. 

Up to this point, we have only considered random removal of links. Another kind of 
removal is targeted removal where the nodes with the largest degree are removed first j^. 
This kind of removal is common in many real world scenarios such as denial of service attacks 
on WWW and delays in airline hubs. 

In scale-free networks, targeted removal of a fraction of q nodes with the largest degree 
can be treated as random removal of q' = g(2-A)/(i-A) ^.j^g network links jsj. After removal, 
the maximum degree is given by K' = mq^^^^~'^\ For A > 3, making the substitutions q q' 
and K ^ K' in Eq. (j4]) we obtain the equation for p^ and the scaling form for Sa (see 
Table [T]). The change to q' and K' reflects the fast collapse of the network and the rapid 
change in the typical network length. The transition line Pc{a) in targeted removal decreases 
significantly more slowly compared to random removal as seen in Fig. [TJ^b). 



In targeted removal for 2 < A < 3, removing even a small fraction of the 
a change in the distance from 2 log log N/ \ log(A — 2) | to log PooN/ log(/t — 1) 
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2^. Thus, 



after percolation Sa can be calculated using the tree approximation which yields 

o log log N n log(K — 1) 

Sa ~ (k- 1) = (log A^) (Scale-free 2 < A < 3, targeted removal). (8) 

In this case, the phase transition to a spanning cluster of order A^ cannot be achieved for 
any finite value of a and p < 1, as seen from Eq. ([8]). Simulation results supporting Eq. ([8]) 
are shown in Fig. [3](d). Comparing random to targeted removal for 2 < A < 3 for LPP yield 
entirely opposite results. In random removal, order A^ nodes are still connected through the 
original paths. On the other hand, in targeted removal for any finite a, the network collapses 
into logarithmically small clusters. 

In summary, our results suggest that the usual percolation theory cannot correctly de- 
scribe connectivity when only a limited set of path lengths are useful. In usual percola- 
tion, order A^ of the network nodes are connected when p > p^. However, in LPP, when 



Pc < p < Pc, only a zero fraction of the network is connected. Therefore, a much smaller 
failure of the network can lead to an effective network breakdown. As an illustration, con- 
sider an ER network with (k) =3, and limit the length between nodes to a = 1.5 times the 
original length. The theory of LPP predicts that the removal of g = 0.31 of the network 
links is enough to break down the network, compared to g = 0.67 in regular percolation. In 
the context of infectious diseases, if the virus typically survives up to 1.5 log steps, our 
theory predicts that the immunization threshold is significantly smaller, 0.31 compared to 
0.67. Due to the above considerations, we expect our results to be important for network 
design, routing protocols and immunization strategies. 
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TABLE I: The functions pc, Sa and (5 for several kinds of network structures under random and 
targeted removal. The scaling of Sa is given for p < pc except for scale-free networks with 2 < A < 3 
under random removal where p is always above pc- 
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FIG. 1: (a) Phase diagram for Erdos-Renyi networks of LPP with respect to parameters a and p, 
demonstrating the linear and power law (fractal) phases for Sa ~ N^- (b) Similar phase diagram 
for scale-free networks with A > 3. The two transition lines represent networks with the same 
K. Note the slow decrease of the transition line for targeted removal compared to the transition 
line for random removal. The region between the two lines has a power law (fractal) phase for 
targeted removal and a linear phase for random removal. In both (a) and (b) the regular percolation 
threshold is fiven by the limit a ^ oo, i.e, Pc = {k)"^ for ER and Pc = {ko — 1)~^ for SF with 
A > 3. 
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FIG. 2: Simulation results (symbols) for Sa vs. N for various network types under random or 
targeted removal and different values of a and p (indicated in plot legends), compared to the 
theoretically predicted power laws (solid lines), with 5 calculated from Table [H Network sizes are 
typically between 1600 and 204800. In all plots the simulation results agree with the theoretical 
predictions, (a) ER networks (random) with (k) = 3 for fixed p = 0.7 and different a values. Inset 
shows the same networks with fixed a = 1.1 and p = 0.5(O)i 0.6(n) and 0.7 ({>). (b) SF networks 
(random) with A = 3.5, m = 2, p = 0.7 and different a values. Inset shows the same networks with 
fixed a = 1.1 and p = 0.5(0), 0.6(n), 0.7(^) and 0.8 (A), (c) SF networks (random), A = 2.2,2.3 
and 2.4, m = 3, p = 0.4, and a = 1. (d) SF networks (targeted) with A = 3.5, m = 3, fixed p = 0.92 
and different a values. Inset shows the same networks with fixed a = 1.2 and p = 0.92(0) and 
0.94 (□). 
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FIG. 3: Simulation results ((a) through (c)) for the scaling of Sa/ic{p)N^) vs. PooN/{c{p)N^) for 
various types of networks, for random and targeted removal for N between 1600 and 25600, and a 
between 1.0 and 4. (a) Erdos-Renyi networks (random) with (k) = 3 and p = 0.5, 0.6 and 0.7. (b) 
Scale-free networks (random) with A = 3.5, m = 2 and p = 0.6, 0.7 and 0.8. (c) Scale-free networks 
(targeted) for A = 3.5, m = 3 and a between 1.01 and 3.0, and targeted removal with p = 0.92 and 
0.94. (d) Simulation results of log Sa vs. log log N and comparison to the theoretical prediction 
(line) for S for SF networks (targeted) for A = 2.3, m = 3, a = 1.5 and p = 0.97. 
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